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Abstract. The main objective of this article is part of a research pro- 
gram to hnk the dynamics of fluid flows with the structure and its tran- 
sitions in the physical spaces. As a prototype of problem and to demon- 
strate the main ideas, we study the two-dimensional Rayleigh-Benard 
convection. The analysis is based on two recently developed nonlinear 
theories: geometric theory for incompressible flows [10] and the bifur- 
cation and stability theory for nonlinear dynamical systems (both finite 
and infinite dimensional) [9J. We have shown in that the Rayleigh- 
Benard problem bifurcates from the basic state to an attractor Ar when 
the Rayleigh number R crosses the first critical Rayleigh number Rc for 
all physically sound boundary conditions, regardless of the multiplicity 
of the eigenvalue Rc for the linear problem. In this article, in addition 
to a classification of the bifurcated attractor Ar, the structure and its 
transitions of the solutions in the physical space is classified, leading to 
the existence and stability of two different flows structures: pure rolls 
and rolls separated by a cross the channel flow. It appears that the struc- 
ture with rolls separated by a cross channel flow has not been carefully 
examined although it has been observed in other physical contexts such 
as the Branstator-Kushnir waves in the atmospheric dynamics [UIZ]- 



The Rayleigh-Benard convection problem was originated in the famous 
experiments conducted by H. Benard in 1900. Benard investigated a fluid, 
with a free surface, heated from below in a dish, and noticed a rather regu- 
lar cellular pattern of hexagonal convection cells. Based on the pioneering 
studies by Lord Rayleigh [I3|, the convection would occur only when the 
non-dimensional parameter, called the Rayleigh number. 



exceeds a certain critical value, where g is the acceleration due to gravity, a 
the coefficient of thermal expansion of the fluid, /3 = \dT/dz\ = (Tg — Ti)/h 
the vertical temperature gradient with Tq the temperature on the lower 
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surface and Ti on the upper surface, h the depth of the layer of the fluid, 
K the thermal diffusivity and the kinematic viscosity. There have been 
intensive studies for this problem; see among others Chandrasekhar [2] and 
Drazin and Reid ^ for linear theories, and Kirchgassner [6], Rabinowitz [12j , 
and Yudovich [16^ I17j. and the references therein for nonlinear theories. 

Recently, the authors have developed a bifurcation theory [9J for nonlinear 
partial differential equations, which has been used to develop a nonlinear 
analysis for the Rayleigh-Benard convections [8|. This bifurcation theory 
is centered at a new notion of bifurcation, called attractor bifurcation for 
nonlinear evolution equations. The main ingredients of the theory include 
a) the attractor bifurcation theory, b) steady state bifurcation for a class of 
nonlinear problems with even order non-degenerate nonlinearities, regardless 
of the multiplicity of the eigenvalues, and c) new strategies for the Lyapunov- 
Schmidt reduction and the center manifold reduction procedures. 

In particular, based on this bifurcation theory, for the Rayleigh-Benard 
convection problem, we have shown [8l |9] that the problem bifurcates from 
the trivial solution to an attractor Ar when the Rayleigh number R crosses 
the first critical Rayleigh number Rc for all physically sound boundary con- 
ditions, regardless of the multiplicity of the eigenvalue Rc for the linear 
problem. 

The main objectives of this article are 1) to classify the solutions in the 
bifurcated attractor Ar, and 2) to study the structure and its transition 
of the solutions of the Benard problem in the physical space. The first ob- 
jective is an important part of the above mentioned new bifurcation and 
stability theory. The second objective is part of a program recently initi- 
ated by the authors to develop a geometric theory of two-dimensional in- 
compressible fluid flows in the physical spaces; see [10]. This program of 
study consists of research in two directions: 1) the study of the structure 
and its transitions/evolutions of divergence- free vector fields (kinematics), 
and 2) the study of the structure and its transitions of velocity fields for 2-D 
incompressible fluid flows governed by the Navier-Stokes equations or the 
Euler equations. The study in this article is in the second direction, linking 
kinematics to dynamics. 

To demonstrate ideas, in this article, we only consider the two-dimensional 
Benard convection problem. The three-dimensional case is technically more 
involved, and shall be reported elsewhere. From the physical point of view, 
two-dimensional Boussinesq equations can be considered as idealized models 
for many physical phenomena, including 1) the Walker circulation over the 
tropics [H], which has the same topological structure as the cells given in 
Figure 14.31 in Theorem 14.21 and 2) the Branstator-Kushnir waves in the 
atmospheric dynamics [H [7|, which have similar topological structure as 
given in Figure ^^2] in Theorem 14.11 

We end this introduction with a few remarks. First, the main idea of the 
study is to explicitly reduce the bifurcation problem to the center manifold, 
together with an attractor bifurcation theorem and structural stability 
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theorems for 2D incompressible flows. The types of solutions in this S 
attractor depend on the boundary conditions. With the periodic boundary 
condition (j2.7p in the xi direction in this article, the bifurcated attractor 
consists of only steady states. When the boundary conditions for the velocity 
field are free slip boundary conditions and Q = (0, L)^ x (0, 1) with < 
< (2 — \/2)/{\/^— 1), using the same method proved in this article, we 
can prove that the bifurcated attractor is still an S^, consisting of exactly 
eight singular steady states (with four saddles and four minimal attractors) 
and eight heteroclinic orbits connecting these steady states. The bifurcated 
attractor and its detailed classification provide a global dynamic transitions 
in both the physical and phase spaces. 

Second, the method and ideas presented in this article are crucial to obtain 
these results, which can not be obtained using only the classical bifurcation 
theories. For the case studied in this article, the classical bifurcation theory 
with symmetry arguments implies that the bifurcation attractor in the main 
theorems. Theorems 4.1 and 4.2, contain a circle of steady states. We need, 
however, the new bifurcated theory to prove in particular that the bifurcated 
attractors are exactly an S^. Furthermore, for general boundary conditions 
such as the free-slip boundary conditions mentioned above, no symmetry 
can be used, and the classical amplitude equation methods fails to derive 
the dynamics. 

Third, the newly developed geometric theory for incompressible flows is 
crucial for the structure and its stability of the solutions in the physical 
spaces obtained in the main theorems. 

Fourth, it appears that the structure with rolls separated by a cross chan- 
nel flow has not been carefully examined although it has been observed in 
other physical contexts such as the Branstator-Kushnir waves in the atmo- 
spheric dynamics [U [7] . 

Finally, as mentioned before, this article is part of a research program 
initiated in the mid 90s to make connections between the dynamics and the 
structure in the physical spaces. 

This article is organized as follows. First the functional setting and the 
attractor bifurcation theorem for the Benard convection problem obtained in 
[8] are introduced in Section 2. Section 3 recapitulates 1) the approximation 
of the center manifold function, 2) attractor bifurcation theorem, and 3) 
structural stability theorems for incompressible flows. The main theorems 
of this article are stated in Section 4, and proved in Section 5. 



2. Benard Problem 

2.1. Boussinesq equations. The Benard problem can be modeled by the 
Boussinesq equations. In this paper, we consider the Benard problem in 
a two-dimensional (2D) domain M} x (0,/i) C (/i > 0). The Boussi- 
nesq equations, which govern the motion and states of the fluid flow, are 
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as follows; see among others Rayleigh |13j . Drazin and Reid ^ and Clian- 
drasekhar [2]: 

du — 

(2.1) _ + (n . V)n - uAu + p^^Vp = -gk[l - a{T - Tq)], 

dT 

(2.2) — + (n-V)r-KAr = 0, 

(2.3) div -u = 0, 

where u,K,a,g are constants defined as in (jl.ip . u = {ui,U2) the velocity- 
field, p the pressure function, T the temperature function, Tq and Ti con- 
stants representing the lower and upper surface temperatures at X2 = 0, /i, 
and k = (0, 1) the unit vector in xs-direction. 
To make the equations non-dimensional, let 

X = hx', 

t = h^t'/K, 

U = KU /h, 

T = ph{T'/VR) + fo- Phx'2, 

p = poK^p'/h^ +po - gpo{hx'2 + (4)^/2), 

Pj. = V I K. 

Here the Rayleigh number R is defined by (II. ip . and = v j k is the Prandtl 
number. 

Omitting the primes, the equations (I2.ip - (j2.3p can be rewritten as follows 



— + (n • V)u + Vp 
ot 



Au - VRTk = 0, 



dT 

(2.5) 'dt^^''' ~ ^""^ - AT = 0, 

(2.6) div u = 0. 

The non-dimensional domain is O = M"*^ x (0,1) C M?. We consider 
periodic boundary condition in the xi-direction 

(2.7) {u,T){xi,X2,t) = {u,T){xi + kL,X2,t) VA; G Z. 

At the top and bottom boundary (x2 = 0,1), different combinations of 
top and bottom boundary conditions are normally used in different physical 
setting such as rigid-rigid, rigid-free, free-rigid, and free-free. For instance, 
we have 

DiRICHLET BOUNDARY CONDITION (rIGID-RIGID) : 

(2.8) T = 0, u = at 2:2 = 0, 1. 
Free-free boundary condition: 

(2.9) T = 0, U2 = ^ = at X2 = 0, 1. 

0x2 
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Free-rigid boundary condition: 

' T = 0, u = at = 0, 

(2.10) < dui 

r = 0, U2 = 0, TT-^ = at X2 = 1. 

The initial value conditions are given by 

(2.11) {u, T) = (uo, To) at t = 0. 

2.2. Functional setting. For simplicity, we proceed in this article with the 
set of boundary conditions given by (|2.7|) and (|2.8|) . and similar results hold 
true as well for other combinations of boundary conditions. 

Let 

H = {{u,T) e L^inf \ divii = 0, ^^21x2=0,1 = 0, ui is periodic in xi-direction}, 
V = {{u,T)e Hl(Slf I divu = 0, (n,T) is periodic in xi direction}, 
Hi = VnH'^{Qf. 

Let G:Hi^H, and Lx = -A + Bx : Hi ^ H he defined by 

Gi^P) = i-P[{u.V)u],-iu■V)T), 
All) = {-P{Au),-AT), 
Bx^ = \{P{Tk),U2), 

for any if: = {u,T) G Hi. Here A = VR, and P the Leray projection to 
fields. 

Then the Boussinesq equations (|2.4p ~ (j2.8p can be rewritten in the follow- 
ing operator form 

(2.12) ^ = L;,V + G(V'), i^ = {u,T). 

2.3. Attractor bifurcation of the Benard problem. Let {Sx{t)}t>o be 
an operator semi-group generated by the equation (|2.12p . Then the solution 
of (|2.12p can be expressed as 

iP{t,i>o) = Sx{t)i'o, t>0. 

Definition 2.1. A set T, C H is called an invariant set of ^2.1^) if S{t)T, = 
T, for any t > 0. An invariant set T, C H of 12.1^) is called an attractor 
if T, is compact, and there exists a neighborhood U d H of such that for 
any ^po £ U we have 

lim distHi^itjtpo),^,) = 0. 

>oo 

Definition 2.2. (1) We say that the equation \2.12\) bifurcates from 
{ip,X) = (0, Ao) to invariant sets ^Ix, if there exists a sequence of 
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invariant sets {^x„} of 112. 12\) such that ^ ilA„ o-nd 
lim A„ = Ao, 

n— +00 

lim max |x| = 0. 

(2) // the invariant sets are attractors of i2.12\) . then the bifurcation 
is called attractor bifurcation. 

We are now in position to state the attractor bifurcation theorem for the 
Benard problem (j2.4p - (j2.8p . The linearized equations of (j2.4p - (j2.6p are given 
by 

' - Au + Vp- VRTk = 0, 

(2.13) I -AT-^U2 = 0, 

, divM = 0, 

where R is the Rayleigh number. These equations are supplemented with 
the same boundary conditions (j2.7p and (j2.8p as the nonlinear Boussinesq 
system. This eigenvalue problem for the Rayleigh number R is symmetric. 
Hence, we know that all eigenvalues Rk with multiplicities of (|2.13p with 
(j2.7p and (|2.8p are real numbers, and 

< Ri < • • • < Rjg < Rk+i < • • • . 

The first eigenvalue -Ri is a function of the period L. The critical Rayleigh 
number Rc is given by 

(2.14) Rc = mmRi{L). 

Let the multiplicity of Rc be mi = m {m = even), and the first eigenspace 
be denoted by Eq. Then we have the following attractor bifurcation theorem. 

Theorem 2.3. [SI [9] For the Benard problem ^2.4\i2.^) . the following asser- 
tions hold true. 

(1) When R < R^, the steady state {u,T) = is a globally asymptotically 
stable in H . 

(2) The equations bifurcate from {{u,T),R) = (0,i?c) to attractors 
for R > Rc, with m — 1 < dimE^^ < m, and S/j is an (m — 1) 
dimensional homological sphere, i.e. has the same homology as 

(3) For any (u, T) £ S/j, the velocity field u can be expressed as 

m / ™ \ 

■u = ^ OfcCfc + o ^ afeefc , 

k=l \k=l / 

where e^ are the velocity fields of the first eigenvectors in Eq. 

(4) For any open bounded neighborhood U C H of {u, T) = 0, the at- 
tractor E/j attracts U \ T in H, where T is the stable manifold of 
{u,T) = with co-dimension m in H. 
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Remark 2.4. Results similar to this attractor bifurcation theorem hold 
true as well for the 3D Benard problems; see [H Ej- 

3. Preliminaries 

3.1. Center manifold functions. To study the structure of the bifurcated 
attractors of (j2.4ll2.8p , it is necessary to consider the reduction of nonlinear 
evolution equations to center manifolds. To this end, we introduce in this 
section a method to derive a first order approximation of the central manifold 
functions, which was introduced and used in ^9j. 

Let H and Hi be two Hilbert spaces, and let Hi ^ H he a dense and 
compact inclusion. We consider the following nonlinear evolution equation 

Lxu + G{u,X), 



(3.1) 



dt 

u(0) = Uq, 



where u : [0, oo) H is the unknown function, A G M is the system parame- 
ter, and L\ : Hi H are parameterized linear completely continuous fields 
depending continuously on A G M^, which satisfy 



(3.2) 



Lx = —A + Bx is a sectorial operator, 

A : Hi H a linear homeomorphism, 

Bx ■ Hi H parameterized linear compact operators. 



It is easy to see (SKUj that Lx generates an analytic semi-group {e~ ^}t>o. 
Then we can define fractional power operators for any < a < 1 with 
domain Ha = D{L'^) such that Ha^ C Ha2 if ai > 02, and Hq = H. 

Furthermore, we assume that the nonlinear term G(-,A) : H^ H, for 
some < a < 1, is a family of parameterized bounded operators (r > 1) 
continuously depending on the parameter A G M^, such that 

(3.3) G{u,X)=o{\\u\\hJ, VAgM\ 

In this paper, we are interested in the case where Lx = —A + Bx are 
sectorial operators such that there exist an eigenvalue sequence {pk} C 
and an eigenvector sequence {e^, hk} C Hi of A: 

' Azk = PkZk, Zk = ek + ihk, 

(3.4) < Repk 00 {k ^ 00), 

^ \lm.pk/iaRepk)\ < c, 

for some a, c > 0, such that {e^, hk} is a basis of H. 

Condition (j3.4p implies that A is a sectorial operator. For the operator 
Bx ■ Hi H, we also assume that there is a constant < 6 < 1 such that 

(3.5) Bx:Hg — > bounded, V A G 



Under conditions (|3.4p and (j3.5p . the operator = — ^ + -Ba is a sectorial 
operator. 
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Let Hi and H be decomposed into 
(3.6) 



Hi = El e El 



\h = eI® E^ 

for A near Aq G M"*^, where E^, E2 are invariant subspaces of Lx, such that 

dimii^^ < 00, 
E^ = El 

E2 = closure of E2 in H. 

In addition, Lx can be decomposed into Lx = Ci ® C2 such that for any A 
near Aq, 



(3.7) 



where all eigenvalues of £2 possess negative real parts, and the eigenvalues 
of C\ possess nonnegative real parts at A = Aq- 

Thus, for A near Aq, equation (jS.ip can be written as 

( dx 



(3.8) 



£iX + Gi(x,y,A), 
^ = £^y + G2(x,y,A), 



where u = x + y G Hi, x G E^ y G E^ Gi{x,y,X) = PiG{u,X), and 
Pi : H ^ Ef are canonical projections. Furthermore, let 

eIo) = closure of E2 in Ha, 



with a < 1 given by (j3.3p . 

The following center manifold theorem is classical; see O [15] . 

Theorem 3.1. Assume h3. 3|) ^( |X7| ) hold true. Then there exists a neigh- 
borhood of Aq given 6y | A — Ao| < <5 for some 6 > Q, a neighborhood Bx G E^ 
of x = 0, and a function $(-,A) : Bx — > E^O) depending continuously 
on \, such that 

(1) $(0,A) = 0, $^(0,A) = 0, 

(2) the set 

Mx = [{x,y)&H\x&Bx, y = $(x. A) G E^e)} , 

called the center manifolds, are locally invariant for \3.1\) . i.e. for 
each uq G Mx 

ux{t,uo)€Mx, VO<t<t(no) 

for some t{uo) > 0, where ux(t,uo) is the solution of liS. 1]) : 
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(3) if {xx{t),yx{t)) is a solution of 113. 8\) . then there is a I3x > and 
kx > with kx depending on (xx{0),yxiO)) such that 

\\yx{t)-'^ixx{t),X)\\H <kxe-f''. 

Also, it is classical that to bifurcation problem of ()3.ip is reduced to that 
for the following finite dimensional system: 

(3.9) ^=Cix + giix,^x{x),X), 
for X e BxC E^. 

Now we give a formula to calculate the center manifold functions. Let the 
nonlinear operator G be given by 

(3.10) G(n,A) = Gfc(n,A)+o(|u|*^), 
for k>2, where Gk{u,X) is a /^-multilinear operator: 

Gk : Hi X ■ ■ ■ X Hi ^ H, 
Gfc(M, A) = Gk{u,--- ,u,X). 
The following theorem was proved in [9]. 

Theorem 3.2. Under the conditions ( f3.3)) -( 3. 7) and /i3.1(]\) . the center man- 
ifold function $(x,A) can be expressed as 

(3.11) ^x, A) = {-C^r'P2Gk{x, A) + 0(|Re/3(A)| • ||xf ) + o(||xf ), 



where C2 is given by ( [g. 7[ ), P2 : H ^ E2 the canonical projection, x G E^, 
and /5(A) = (/3i(A),--- ,/3m(A)) the eigenvalues o/£^. 

Remark 3.3. Consider the case where Lx ■ Hi H is symmetric. Then 
the eigenvalues are real, and the eigenvectors form an orthogonal basis of 
H. Therefore, we have 

u = x + y £ El eE2, 



X = ^XjCj G El, 



i=l 

00 

y = ^ X'lei G E2. 

i=m+l 

Then near A = Aq, the formula (j3.1ip can be expressed as follows. 

00 

(3.12) ^{x,X)= + 0(|Re/3(A)| • ||xf ) + o(||a;f ) 

j=m+l 

where 



/3,(A) ^ -h-h''h---^h, 
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In many applications, the coefficients Oj^.-.j^, can be computed, and the first 
m eigenvalues /3i(A), • • • satisfy 



|Re/3(Ao) 



\ 



^(Re/3,(Ao))2 = 0. 



Hence (j3.12p gives an explicit formula for the first approximation of the 
center manifold functions. 



3.2. S'^-attractor bifurcation. In this section, we study the structure of 
the bifurcated attr actor of (j3.9p when m = 2. Namely, we consider a two- 
dimensional system as follows: 



(3.13) 



^ = /3(A)x-5(x,A), xGM^ 



Here /3(A) is a continuous function of A satisfying 



(3.14) 



mi 



( <0 

= 
> 



if A < Ao, 
if A = Ao, 
if A > Ao, 



and 



(3.15) 



' g(x,\) = gk(x,X)+o{\x\^), 

gk{-, A) is a k-multilinear field, 
.Ci|x|^+^ < (5fc(x,A),x), 



for some integer k = 2m + 1 > 3, and some constants < Ci. 

The following theorem was proved in which shows that under condi- 
tions (|3.14p and (|3.15p . the system (|3.13p bifurcates to an 5^-attractor. 
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Figure 3.1. has AN + n {N = 1 and n = 2 shown here) 
singular points, where pi and p4 are saddles, ps and pq are 
nodes, and p2 and p^ are singular points with index zero. 

Theorem 3.4. Let the condition Jg. i^p and ii^.l5\) hold true. Then the 
system \3.1'J\) bifurcates from (x,A) = (0, Aq) to an attractor Tj\, which is 
homeomorphic to S^, for Aq < A < Aq + e for some e > 0. Moreover, one 
and only one of the following is true. 

(1) Sa is a periodic orbit, 

(2) Tix consists of infinite number of singular points, or 

(3) Tix contains at most 2{k + l) = 4(m + l) singular points, consisting of 
2N saddles points, 2N stable node points, and n (< A{m + 1) — AN) 
singular points with index zero, as shown in Figure \3.1\ 

3.3. Structural stability theorems. In this subsection, we recall some 
results on structural stability for 2D divergence-free vector fields developed 
in [inj, which are crucial to study the asymptotic structure in the physical 
space of the bifurcated solutions of the Benard problem. 

Let C"'(1^,M2) be the space of all C (r > 1) vector fields on 1] = x 
(0, 1), which are periodic in xi direction with period L, let D'^{U,, M^) be the 
space of all C divergence-free vector fields on 17 = x (0,1), which are 
periodic in xi direction with period L, and with no normal flow condition 
in a;2-direction: 



D'\n, M2) = {ve C""(17, M^) I W2 = at X2 = 0, 1 } . 



Furthermore, we let 



Bl{VL, R2) = {v£ D''{n, R2) I t; = at a;2 = 0, 1 } , 
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Definition 3.5. Two vector fields n, t> G C"^(0,M^) are called topologically 
equivalent if there exists a homeomorphism o/ (/9 : $7 — > fi, which takes the 
orbits of u to orbits of v and preserves their orientation. 

Definition 3.6. Let X = D''{n,m?) or X = B^{n,M.'^). A vector field 
V & X is called structurally stable in X if there exists a neighborhood U C X 
of V such that for any u €z U , u and v are topologically equivalent. 

Let V G D''(0,R2). We recall next some basic facts and definitions on 
divergence-free vector fields. 

(1) A point p £ Q is called a singular point of v if v{p) = 0; a singular 
point p of is called non-degenerate if the Jacobian matrix Dv{p) 
is invertible; v is called regular if all singular points of v are non- 
degenerate. 

(2) An interior non-degenerate singular point of v can be either a center 
or a saddle, and a non-degenerate boundary singularity must be a 
saddle. 

(3) Saddles of v must be connected to saddles. An interior saddle p ^VL 
is called self-connected if p is connected only to itself, i.e., p occurs 
in a graph whose topological form is that of the number 8. 

Let V G Bq{Q,, M^); then we know that each point on X2 = 0, 1 is a singular 
point of V in the usual sense. To study the structure of v, we need to classify 
the boundary points as follows. 



Definition 3.7. Let u G Bl{TM){r > 2). 

dn 



1) A point p G DM is called a d-regular point of u if ^"T*-^^ ^ 0; other- 



wise, p G dM is called a d-singular point of u. 
A d-singular point p G dM of u is called non- degenerate if 



/ d'^Urip) d'^Ur{p)\ 

(3.16) det 



drdn dn^ 



d'^Unjp) d'^Unjp) 
V drdn dn^ ' 



A non- degenerate d-singular point of u is also called a d-saddle point 
ofu. 

o 

(3) u G Bq{TM) (r > 2) is called D-regular if a) u is regular in M, and 
b) all d-singular points of u on dM are non- degenerate. 

The following theorem provides necessary and sufficient conditions for 
structural stability of a divergence-free vector field. 

Theorem 3.8. ^ Let u e Bl{TM){r > 2). Then u is structurally stable 
in Bq{TM) if and only if 

1) u is D-regular; 

2) all interior saddle points ofu are self- connection; and 
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3) each d-saddle point of u on dM is connected to a d-saddle point on 
the same connected component of dM. 

Moreover, the set of all structurally stable vector fields is open and dense in 



Remark 3.9. For vector fields with free-rigid boundary conditions, the con- 
ditions for structural stability differs slightly. More precisely, u £ B'[{TM){r > 
2) is structurally stable in B\{TM) if and only if 

1) all singular of n in ^7 and on 3:2 = 1 are regular, and all (9-singular 
points on X2 = are 9-regular; 

2) all interior saddle points of u are self-connected; and 

3) each saddle of u on 2:2 = 1 is connected to saddles on X2 = 1, and 
each 9-saddle point of u on X2 = is connected to a 5-saddle point 
on X2 = 0. 

Remark 3.10. For vector fields satisfying free-free boundary conditions, 
we set 



Then u G B^{n,R'^) (resp.n e Bl(n,R'^)) is structurally stable in B^{n,R'^) 
(resp. in i?2(J^,M^)) if and only if 

1) u is regular; 

2) all interior saddle points of u are self-connected; and 

3) each boundary saddle of u is connected to boundary saddles on the 
same connected component of dQ (resp. each boundary saddle of u 
is connected to boundary saddles not necessarily on the same con- 
nected component). 

The difference between these two cases is due to the zero-average condition in 
the definition in S3(r2,M^), which implies that B^{Q,R'^) does not contain 
the harmonic field vq = (a, 0) for any constant a ^ 0. Hence, an orbit 
connecting two saddles on different components of the boundary can not be 
broken with a perturbation in i?3(r2,M^) into orbits connecting only saddles 
on the same connected component of the boundary. 

4. Structure of Bifurcated Solutions for the Benard Problem 

In this section, we study the topological structure of the bifurcated at- 
tractor and the asymptotic structure of solutions for the Benard problem. 
It is known that for each type of boundary conditions, there is a minimal 
period Lc satisfying (I2.14p . Hereafter, we always take Lc to be the period 



Bl{TM). 




— — = at X2 = 0, 1 





of ([2TD . 



The main theorem in this article is as follows. 
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Figure 4.1. All points on the bifurcated attr actor J^ji 
are steady state solutions. 



5^ 



Theorem 4.1. For the Benard problem ^2.4^ - [¥7^) . the following assertions 
hold true. 



(1) For R > Rc, the equations bifurcate from the trivial solution {{u,T),R) 
(0, i?c) to an attractor S/j, homeomorphic to , which consists of 
steady state solutions as shown in Figure\4jJ\ where Rc is the critical 
Rayleigh number. 

(2) For any ipQ = (uq, Tq) G H\{T L) E), there exists a time to ^ such 
that for any t > tg, the vector field u{t, tpo) is topologically equivalent 
to the structure as shown in either Figure \4^ a) or Figure \4^ b), 
where V' = {u{t,'^po),T(t,^jJo)) is the solution ip = (^(t, V'o)) ^(^i "^o)) 
of ^2.4\ )- f278\) with initial data -00; T is the stable manifold of the 
trivial solution (u, T) = with co-dimension 2, and 



E = Uu,T) e H 



uidx2 = 
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(b) 

Figure 4.2. Here the horizontal axis is the xi-axis, and the 
vertical axis is the X2-axis. With the Dirichlet boundary 
conditions on 3:2 = 0, 1, the flow is not moving on both the 
top X2 = 1 and the bottom X2 = boundaries. 

The zonally moving meandering flow shown in Figure 14.21 appears often 
in many physical problems such as the Bransdator-Kushnir waves in atmo- 
spheric circulation [H [7] . 

Theorem O is also vahd for the Benard problem (fT4]) - ([2?7l) with (f2rT0]) . 
However the case with the free- free boundary condition is different. More 
precisely, for the free-free boundary condition, it is easy to see that for any 
constant q, the harmonic field Vo = ((cK)0),0) is a solution of (I2.4p - (j2.6p . 
Therefore, we have to consider the problem (I2.4p - (|2.7p with ()2.9p in the 
following function spaces: 

H = {{u,T) € L^{nf I divn = 0, / ndx = 0}, 

Jn 

Hi = {{u,T) £ H n H^{nf satisfies ([MD and ^ } 
Then we have the following theorem. 

Theorem 4.2. For the Benard problem \2. 4^ - ^27^ with boundary condition 
^2. 9^) . the following assertions hold true. 

(1) For R > Rc, the equations bifurcate from the trivial solution {{u,T),R) = 
(0,i?c) to an attractor T,r, homeomorphic to , which consists of 
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steady state solutions, where Rc = 277r^/4 is the critical Rayleigh 
number. ^ 
(2) For any TpQ = {uo,Tq) S H\T, there exists a time to > such that for 
any t > to, the vector field u{t,ipo) is topologically equivalent to the 
structure as shown in Figure where = {u{t,tpo),T{t,tpo)) is 
the solution ^ = {u{t, i'o),T{t, -^o)) of [24\)-[2l\) with / fO) . T is the 
stable manifold of the trivial solution (u, T) = with co-dimension 
2 in H. 




Figure 4.3. Here the horizontal axis is the xi-axis, and the 
vertical axis is the X2-axis. With the free slip boundary con- 
ditions on X2 = 0,1, the flow does move on both the top 
X2 = 1 and the bottom X2 = boundaries. 



5. Proof of Main Theorems 



5.1. Eigenvectors of the linear Boussinesq equations. We shall only 
prove Theorem 14. 1[ The proof of Theorem 14.21 is essentially the same, and 
we omit the details. We proceed by first considering the eigenvalues and 
eigenvectors of the linearized equations of (j2.4p - (j2.6p : 

' Au-Vp + VRTk = (3{R)u, 

(5-1) < AT + VRu2 = P{R)T, 

, divu = 0, 

supplemented with the boundary conditions (j2.7p and ()2.8p . 

For ip = (mi,U2,T) G Hi, we take the separation of variables as follows: 



2A;7rxi , 2A;7r 2kTTXi 2kTTXi , 

sm — - — h {x2), —— cos — - — n(x2), cos — - — 0{x2, 



cos ■ 



L 

2kTTXl 



L 

, 2kTT 

h [X2), —r- sm 



2kTTXi 2k-KXi , 

h{x2),sm — - — e{x2. 



L ' L L ' ^' L 

Then it follows from (jS.ip that (/i, 6) satisfies the following differential equa- 
tions 

J2 \ 2 / ,2 



(5.2) 



\dxl 
_ 



h 



RokB 



dxo 



2 

0-k 



Rauh 



= -P{R)e, 



h, 
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supplemented with the following boundary conditions 

(5.3) = o,h = h' = O atx2 = 0,l, 

where = Ikn/L and L = satisfies (I2.14p . 

The eigenvalue problem (j5.2p with (j5.3p is symmetric, and has a complete 
eigenvalue and eigenvector sequences for given k and R: 

( Pki{R) > l3k2{R) > ■ ■ ■ , 



(5.4) 



Moreover, 



lim (3kj{R) = -oo, 

hkj£H\0,l)nHiiO,l) i = l,2,--- 

ekjeH\o,i)nH^{o,i) j = i,2,--- 



{{hkj,6kj) I J = 1,2,---} 
constitutes an orthogonal basis of iv^(0, 1) x L^(0, 1). 

Thus, we obtain the following complete set of eigenvectors for ()5.ip with 
boundary conditions ()2.7p and ()2.8p : 

(5.5) 

, / . 2k'KXi . , , , Ik-K Ikixxx , , , 2kTxx\ „ , A 

V'fcj = |-Sm ^ /tfc^(X2), — cos ^ %j(X2),COS ^ 6'fcj(X2) I , 

(5.6) 

~ / 2k-KX\ , , , , 2A;7r . 2kTTXi .... 2A;7rxi ^ , A 
V'fcj = I cos — - — h,^j{x2),^--sm — - — hkj{x2),sm — - — 0kj{x2) I . 

where < A; < oo, 1 < j < oo. When A: = 0, we derive from (j5.ip . (j5.5p and 
(15:6]) that 



(5.7) 



r -^o^- = (0,0,sinj7rx2), 
[i^Oj = ih'ojix2),0,0). 

5.2. Singularity Cycle. We shall show that the bifurcated attractor S/j of 
(j2.4p - (j2.8p given in Theorem 12.31 contains a cycle of steady state solutions. 

First, we note that the equations (|2.4p - ()2.8p are invariant under the fol- 
lowing translation: 

ip{xi,X2,t) ^ ip{xi + a,X2,t) Va € M. 

Hence, if ipoix) is a steady state solution, then ^0(2:1 + a,X2) are steady 
state solutions as well. It is easy to see that the set 

S = {'ipo{xi + a, X2) I a e M} 

is a cycle in Hi. Therefore, each steady state of (|2.4|) - (|2.8|) generates a 
cycle of steady state solutions. 
Let 

H' = {(n,T) G H I ni(-xi,X2) = -ni(xi,X2)}, 
H[ = Hin H'. 
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It is easy to check that H' and are invariant spaces for the operator 
L\ + G given by (I2.12P in the sense that 

L\ + G : H'l H', 

where A = On the other hand, it is clear that the sequence {ipkj I ^ = 
0, 1,-- - ,j = 1,2, •••} defined by (15. Sh is a basis of H'. Since the first 
eigenvalue Rc of (j2.13p is simple, the first eigenvalue Pi{Rc) of Lx in H[ is 
also simple, where 

Pi{Rc) = /3u{Rc) = 0, 

and (3ii{Rc) is defined by ()5.4p . Hence by the classical Krasnoselskii bifur- 
cation theorem, we know that the operator Lx + G bifurcates from (ip, A) = 
(0, \/^) to a singular point in H[. Namely, the Benard problem (j2.4p - (j2.8p 
bifurcates from (ip, R) = (0, Rc) a steady state solution. Therefore, the 
bifurcated attractor T,ji contains at least a cycle of steady state solutions. 



5.3. S'l-attractor: = S^. To prove that = S'^, by Theorem [331 
it suffices to verify that the reduced equations of (j2.4p - ()2.8p to the center 
manifold satisfy conditions (I3.14P and ()3.15p . 

For any -0 = (u, T) G H, we have 

k>0,j>l 

Since L is the minimal period satisfying (12.14p . ipn and ipu are the first 
eigenvectors of ()5.ip . Therefore the reduced equations of ()2.4p - (l2.8p are 
given by 



(5.8) 



( ^ = (3i{R)xn + ^-^(G(0,0),V.n), 



-/3i(i?)yii + ^;^;-^(G(V,V),V'ii). 



Here for Vi = (^t,Ti), 02 = (t^,72), and ^-^ = {w,T^), 

(G(0i,02),03) = - / [{u-Vv)w + {u-VT2)Ti]dx. 
Jn 

Let the center manifold function be denoted by 

(5.9) ^ (^^kjixu,yu)ipkj + ^kj{xii,yii)^pkj^ . 

Note that for any il^i e Hi [i = 1,2,3), 

(G(0i,02),V'2) =0, 

(G(0i,02),V'3) = -(0(^1,03,), 02). 



RAYLEIGH-BENARD CONVECTION 19 

Then hy ip = xn-i/;!! + yiiV'ii + we have 

(5.10) =(G(V^ii,V^ii),Vii)y?i 

- {G{ipu,ipn),'>Pn)xnyu 

- (00^11,11^11), <^)xn 

+ + G($,V'ii),V'ii)yii 

+ (G($,$),Vii), 

(5.11) (G(7/;,V),V^ii) =(G(7/;ii,V'ii),V^ii)x?i 

- (G(V'ii,-0ii),'0ii)xiiyii 

- (G(V^ii,V^ii),$)yii 

+ (^(1/^11,$) + G($,i/'ii),V^ii)xii 

It is easy to check that for / 0, 2, 

' (G'(V'll,^ll),'0A:j) = 0, 

(G(V'ii,V'ii),V'fcj) = 0, 
(G(V'fcj,^ii),'0ii) = 0, 
^ (G(V'fcj,'0ii),V'ii) = 0, 
and for any k > 0, j > 1, we have 

' (G(V'ii,V'ii),V'fcj) = 0, 

(G(V'ii,V'ii),V'fcj) = 0, 

(G(V'fcj,V'ii),'0ii) = 0, 

, (G(V'fcj,^ii),V'ii) = 0. 

By dSSt), (f5l^ and (f533|l . the equahties (f530]l and (fETTD can be rewritten 
as 

(5.14) 

00 

(Giip, V;), V'li) = - J]][(G(^ii, V'li), V'oj)'^oj + (G'(V'ii, -011), ^2j)^2i]xii 



(5.12) 



(5.13) 



i=i 



^[(0(^^11, ^11), V'oj) + (G(V'oj,'0ii),V'ii)yii^. 



00 



-Yl[{Gi'lpu,iJu),1p2j) + (G(^2i,V'll),V'll)yil^2i 

+ (G($,^'),Vii), 
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(5.15) 

oo 

{G{'tp,^),^ii) = - ^[(G(^ii,V^ii),V'oi)$oj + (G(V^ii,V^ii),V'2j)^2j]yii 



oo 



^[(G(V'ii,'i/'ii),V'oj) + (G(V'oj,V'ii),V'ii)a;ii^oj 



oo 

~ ^[{G{'<Pii,'ifii),'ip2j) + (G(V'2j,V'ii),V'ii)a;ii52j 

i=i 

+ (G($,$),V^ii). 

Since the center manifold functions contains only higher order terms 
Hxiuyii) = 0{\xnf,\ynf), 

we derive that 



(5.16) 



r (G($,$),Vii) = o(|xii|Myn|'), 
\(G($,$),V^n) = o(|xn|Myn|')- 



Then direct calculation yields that 

(5.17) {G{^2j,M,M = - iG{^2j,^u),i^u) 

=ll [-h[Ah'2jh'u + 2h2jh'l,) 
+ hn{h'2jhn + 2h2jh[i) 

+ eiMjen + 2h2je'u)]dx2 

= 2 d^^~^2j(/i'ii)' + h2jhi^ + h2j9li)dx2 
=0. 



It is clear that for any A; > and j >1, 



(5.18) 



UkjfH = UkjfH. 
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Hence the reduced equations (|5.8p can be expressed as follows: 

(5.19) ^ =/3i(i?)xn - — ^ V[(G(Vii,Vii),V'o,)$OiXn 

\mi\\H~[ 

+ (G(V'll,V'll),'02j)^'2jXll + (G(V'll, V'll),'02j)yil^2j 

+ {G{4)u,Tpii),ipoj)yn<^oj + (G(v5oi, V'li), V'ii)yii^'oj] 
+ o(|xii|Myii|3), 

(5.20) ^ =Pi{R)yii - — ^ V[(G(V'ii,^ii),V'o,)%yii 

+ (G(V'll,'01l), V'2j)^2jyil + iG{lpu,i'u),'^2j)xu^2j 

+ (G(V'ii,'0ii),^oj)a;ii^oj + (G(^oi,V'ii),'0ii)a;ii^oi] 
+ o(|xn|Myii|=^). 
By Theorem 13.21 and (|3.12|) , we have 

^oi = 11, 112^ [{G{'ilJu,ipii),ipoj)xn + {G{i)u,ipn),ipoj)yn] 

WVOjllHlJOj 

+ o{xi, + yl) + 0{Pi{R){xl,+yf,)), 

= 11 , 112 ,Q [{G{ipn,i^n),i^2j)xu + (GCV'ii, V'li), '02j)y?i] 

+ o(x?i + yj,) + 0{(3,{R){xl + y?J), 
1 

+ oixl, + yl,) + OmR){xl+yl,)), 
= „~ „■! ^ [(G'(V'ii.'0ii) + (G(V^ii,V'ii),V'2i)2;iiyii] 

+ o(x?i + y?i) + 0(/3i(ii)(x?i+y?i)), 
=o(x2, + VA;/0,2, 

=o(x?i + y?i) VA;/0,2, 

where PkjiR) are as in (15. 4p . 

Also by direct computation, we obtain that 

(G(V'ii,'0ii),'0oi) = (G(V'ii,^ii),'0oi), 
(G(V'ii,V'ii),V'2i) = -{G{i)n,i'ii),i^2j), 
(G(V'ii,V'ii) + G('^ii,'0ii),V'oj) = 0, 

(G(V'll,V'll),^2i) = (G(V'll, V'll),V'2i) = (G(V'll,'01l),V'2j)- 



= „T i|2 ^ [(^'(V'li^V'ii) + (G(V'ii,V'ii),V'oi)a;iiyii] 
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Hence we have 

(5.21) $0, , ~L {G{^Pn,M,^oj){x'n+yl) 

+ o{xl,+yl,)+OmR)ixl,+yl)), 

(5.22) ^^.=^^(^(^n,^n),^2.)(x?,-.L) 

+ o{xl,+yl,)+0{(3i{R){xl, + yl,)), 

(5.23) =o(x?i + + 0(/3i(i?)(x?i + y?i)), 
~ —2 

(5.24) $2i = 11 , 112 ^ (g(^ii,-0ii),^2j)xiij/ii 

+ o(x?i+y?i)+0(/3i(i?)(x?i+y?i)). 
Inserting (fOB - dOi]) into ([5T9]) and ([OO]) . we have 

(5.25) ^ =/3i(i?)a;ii - axuixj^ + yfi) 

+ o(x?i + y?i)+0(/3i(i?)(x?i+y?i)), 

yii - aynixli + yfi) 
+ oixl + yl) + OiPiiR){xl,+yl)), 



(5.26) ^ =/3i(^)yii - ayii(x?i + yj,) 



where 



a 



^ oo 



(G(V;ii,V;ii),^o,)^ ^ {G{i^u,M,i^2jy 



We know that 



(3i{R) = Pn{R) > Pkj{R) "^(kj) + (1,1), 
/3ii(i2c) = 0. 

Hence near R = Rc, 

Pk,{R)<0 V(A:,j) / (1,1). 

Consequently, q > and (^^TIB and (fOUj) satisfy ([XT^ . 
Also, we know that [81 Ej 



/3i(i?)<( 
Therefore (|3.15p holds true. 



r < if R<Rc, 

= iiR = Rc, 

> a R> Rr. 
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5.4. Asymptotic structure of solutions. By [4J, we know that for any 
initial value ipo = {uo,To) G H, there is a time r > such that the solution 
tp = (tt(t, V'o)) V'o)) is C°° for t > T, and is uniformly bounded in C"*- 
norm for any given r > 1. Hence, by Theorem 12. 3^ we have 



(5.27) 



lim min ||V'(t,'0o) - 4>\\c- 



0. 



We infer then from (j5.25p and (j5.26p that for any steady state solution 
$ = (e, T) G Tjr of (j2.4p - (l2.8p . the vector field e = (ei, 62) can be expressed 
as 

27r 

ei = rcos -j^ixi + 6*) /iii(x2) + yn, 
62 = ^rsin^(xi + 6) hu{x2) + V2(xii, yn, 




Vi{xn,yn,l3i) = o{yJ [3i{R)) 



for i = 1,2. 



On the other hand, it is known that the first eigenfunction hii{x2) of (j5.2p 
and (|5.3p at R = Rc is given by 

(5.30) hii{x2) ^ cos ao(x2 ~ 2^ ~ ^-^^ '^'-'^^ 01(^2 "2-^ a2(2;2 "2-^ 
+ 0.1sinhai(x2 - -)sina2(x2 - -), 

where ao — 3.97, ai ~ 5.2 and a2 ~ 2.1. This function is schematically 
given by Figure 15. It see [2j . 




Figure 5.1. 
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Now we show that the vector field 
(5.31) eo ■■ 



211X1 , , , ^ 27rr . lirxi 
r cos — —hii[X2), —;— sm — —hii{X2) 



L 



L 



L 



is D-regular in $7 = x (0, 1). 
To this end, by (j5.30p we see that 

^'/i(x2)/0 at X2 = 0,^,1, 

Hence 

det DeQ{xi,X2) / 0, 
for any (xi, X2) = (A;L/2, 1/2) with /c = 1, 2, • • • , and 
/ d^el d\l\ 

27rr 



det 



dxidx2 dxo 



— sm — ^ /ill (2:2) T 0, 



\ 9X1 0X2 5x2 / 

for any (xi,X2) = ((2A; + l)L/4, 0) or (xi,X2) = {{2k + l)L /A,l) with k = 
1,2,- •• . Therefore, the vector field (j5.3ip is D-regular, and consequently, 
the vector fields e in (j5.28p are D-regular for any Rc < R < Rc + £ ^oi some 
e > small. 

Next we show that the following subspace of H 

E = { {u,T) e Hi \ [ [ uidx = 0} 
Jo Jo 

is invariant for ()2.4p - ()2.8p . In fact, we can verify that 

(5.32) [ [ P[{u ■ V)u] ■ idx = yueHi, 

Jo Jo 

where i = (0, 1)* is the unit vector in the xi-direction, and P the Leray 
projection. Indeed, by the Helmholtz decomposition, we have 

[{u-V)u]-i = P[{u-V)u]-i + -^^, 



for some c/) G H^{Q). Hence, 



^0 



P[{u ■ V)u] ■ idx 



JO 



{u ■ 'V)uidx = 0. 



The invariance of E for ()2.4p - ()2.8p implies that for the vector field e given 
in (|5.28p . we have 



(5.33) 



Jo 



Vidx2 = 0. 
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Hence in the Fourier expansion of e in ()5.5p and (15. 6p . the coefficients of 
■0Oj are zero. By the connection lemma in [10], it follows from (j5.33p that 
the vector field e = (ei, 62) of (|5.28p is topologically equivalent to the vector 
field Co given by (j5.3ip . which has the topological structure as shown in 
Figure liTSl 

For any initial value ipo = {uq,Tq) £ H \ E, 

k 

$0 G E, 

V^Ofc = (sin7ra;2, 0, 0), k = 2m + 1, m = 0, 1, • • • . 
By (|5.32p . we see that for the operator G in (j2.12p . 

(G(V'),V^ofc) = 0, 

which implies that the solution of p.4| )-( |2T8]l with ()2.1ip has the following 
form 

(5.34) V(i,V'o) = 5^afce*^«'=V^ofc + $(t, V'o), 

k 

where ^ £ E, and /3ofe < near R = Rc- 
Let 

K = min{ k \ aj. ^ Q and k is odd}, 

and let ^{t.ipo) = {u{t,ijo),T{t,iJo)) be the solution of jS^ll-dlS]) given by 
(j5.34p . Then by (j5.27p . the vector field n(t, V'o) is topologically equivalent 
to the following vector field for any t > to with to > sufficiently large 

(5.35) n = e + (a^e"*** sinA;7rx2,0), 

where e is as in (j5.28p . 

Since the vector field e is L>-regular and topologically equivalent to the 
vector field as shown in Figure 14.31 Then using the method for breaking 
saddle connections in [lOJ, it is easy to show that the vector field u given by 
(j5.35p is topologically equivalent to the structure as shown in Figure I4.2l fa) 
if Ofc < 0, and to the structure as shown in Figure S^l^b) if a/. > 0, for any 
t > to sufficiently large. 

Thus the proof of Theorem 14.11 is complete. 
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